Resonant nature of phonon-induced damping of Rabi oscillations in quantum dots 
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Optically controlled coherent dynamics of charge (excitonic) degrees of freedom in a semiconductor 
quantum dot under influence of lattice dynamics (phonons) is discussed theoretically. We show that 
the dynamics of the lattice response in the strongly non-linear regime is governed by a semiclassical 
resonance between the phonon modes and the optically driven dynamics. We stress the importance 
of the stability of intermediate states for the truly coherent control. 



One of the most challenging proposals exploiting the 
atomic-like properties of quantum dots (QDs)i*£ is to im- 
plement the quantum information processing schemes' 5 ' 4 
using confined carrier states in these artificial systems. 
One of the proposed solutions is to use the charge de- 
grees of freedom which may be efficiently controlled 
optically even on sub-picosecond scale^SiS which may 
look promising in view of the measured 1 ns lifetime 
of excitons in QDs 910 . The first steps towards this 
goal include experimental proofs of controlled coher- 
ent dynamics in these structures^ 1 -, observation of Rabi 
oscillations^i&iirfiiiiiii, demonstration of entanglement 
between states of interacting dotsii and performing a 
quantum logic gate based on a bi-exciton system in a 
QDi£. 

However, the coherent dynamics of quantum confined 
carrier states is very sensitive to interaction with the 
macroscopic number of degrees of freedom of the out- 
side world. If the coherent dynamics involves an ex- 
cited state of the confined excitonAiii&ii the coherence 
time is limited by relaxation to the ground state (usu- 
ally ~ 30 ps). The measured signal may further be af- 
fected by the particularities of the detection technique^. 
Another approach consists in measuring the dot occu- 
pation upon driving by pulses with constant length and 
varying amplitudesi 2 * 1 ^* 1 ^. So far, it has always turned 
out that such pulse-area-dependent Rabi oscillations de- 
viate from the ideal ones, the discrepancy being larger for 
stronger pulses. In principle, this might also be explained 
by experimental conditions or environmental perturba- 
tion: scattering by weakly localized excitons around an 
interface fluctuation QD (further confirmed by increas- 
ing decay for stronger pulses)^, tunneling to leads in the 
photodiode structure (on ~ 10 ps timescale)^, or dipole 
moment distribution in the QD ensemble 9 . 

One might believe that all the perturbation comes from 
sources that may be removed or minimized by technology 
improvement and by optimizing the experimental condi- 
tions and hence produce no fundamental obstacle to arbi- 
trarily perfect quantum control over the excitonic states. 
However, in every case the QDs are inherently coupled 
to the surrounding crystal lattice. The perturbing effect 
of lattice modes (phonons) has been observed experimen- 
tally as a fast (~ 1 ps) partial decay of optical polariza- 
tion induced by an ultra-fast laser pulse 9 -. The theoretical 
analysisiiiSSiSi shows that the experimentally observed 
effects may be quantitatively accounted for by invoking 
the carrier-phonon coupling. Indeed, under suitable ex- 
citation conditions in a high-quality sample the carrier 



density may be kept low enough to eliminate Coulomb 
scattering effects while the coupling to the electromag- 
netic modes is manifested by the exciton radiative decay 
on timescales many order of magnitude longer than these 
relevant here. The theory shows that the coherence decay 
should be viewed as a trace of coherent lattice dynam- 
ics (due to lattice inertia), persisting even at zero tem- 
perature when the lattice is initially at its ground state 
and cannot perturb the stationary ground state of the 
carrier subsystem. Due to strong reservoir memory on 
timescales relevant for these processes, they cannot be 
fully understood within the Markovian approximations. 
In particular, the idea of a "decoherence time", with 
which the control dynamics competes, is misleading 2 ^. 

The recent theoretical study2£ on optical Rabi flopping 
of excitons in QDs driven by finite-length optical pulses 
shows that exponential damping models fail to correctly 
describe the system kinetics. The appropriate descrip- 
tion yields much less damping, especially for long pulses. 
It turns out that the lowest "quality" of the Rabi oscil- 
lation is obtained for pulse durations of a few ps, while 
for longer durations the damping is again decreased. 

In this paper we propose both qualitative and quantita- 
tive explanation of the mechanism leading to the phonon- 
induced damping reported in the theoretical and exper- 
imental studies. We show that the carrier-phonon inter- 
action responsible for the damping of the oscillations has 
a resonant character: While in the linear limit the system 
response depends only on the spectral decomposition of 
the pulse, the situation is different when a strong pulse 
induces an oscillating charge distribution in the system. 
In a semi-classical picture, this would act as a driving 
force for the lattice dynamics. If the induced carrier dy- 
namics is much faster than phonon oscillations the lattice 
has no time to react until the optical excitation is done. 
The subsequent dynamics will lead to exciton dressing, 
accompanied by emission of phonon packets, and will 
partly destroy coherence of superposition state a 9 i 19 i 20 i 21 
but cannot change the exciton occupation number. In 
the opposite limit, the carrier dynamics is slow enough 
for the lattice to follow adiabatically. The optical excita- 
tion may then be stopped at any stage without any lattice 
relaxation incurred, hence with no coherence loss. The 
intermediate case corresponds to modifying the charge 
distribution in the QD with frequencies resonant with 
the lattice modes which leads to increased interaction 
with phonons and to decrease of the carrier coherence 
(see Ref. 0for a simple, single- mode model). 

The above ideas are supported by the formal analy- 
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sis in the following. First, we describe the theoretical 
method used; then the description is applied to a specific 
self- assembled InAs/GaAs quantum dot to get quantita- 
tive estimations; finally, we conclude with some summary 
and possible consequences of the results. 

Following the earlier studiesiS* 2 ^' 2 ^ successfully ac- 
counting for the observed carrier dynamics on picosecond 
timescales, we describe the system by the RWA Hamil- 
tonian (in the rotating frame) 

H = H x + H ph + H int . 

Here Hx = f{t){a + a') describes the interaction of the 
ground-state exciton with the resonant laser pulse mod- 
ulated by the real envelope function f(t) (a, a' are exci- 
tonic operators), i/ p h = Sk^fe&fe&fc is the lattice Hamil- 
tonian (b k , b k refer to the phonon mode k whose energy 
is Wfe), and 



H int =a^aY J (fkbl + & k 



is the interaction Hamiltonian (f k are the coupling con- 
stants for the individual phonon modes) . Due to momen- 
tum conservation, only long wavelength phonon modes 
(compared to dot size) are effectively coupled 2 ^* 2 ^. For 
the timescales relevant here, the damping effect is caused 
mostly by acoustic phononsS&Si. In the calculations, we 
include deformation potential coupling to longitudinal 
acoustical (LA) phonons, which describes the modifica- 
tion of the band structure due to lattice compression and 
is different for electrons and holes. Since no cancellation 
is involved, its effect depends only weakly on the exact 
wavefunction geometry. We assume that transitions to 
higher states, including possible inter-valley transitions, 
are well separated in energy, compared to the spectral 
characteristics of the dynamics considered here and may 
be neglected. Also bi-exciton transitions are assumed to 
be eliminated by a suitable polarization choice. Piezo- 
electric coupling to acoustic phonons has been shown to 
be negligible if the electron and hole function overlap^ 
and is neglected here (but might contribute if the wave- 
function geometry is different). We have verified that 
the Frohlich coupling to the longitudinal optical (LO) 
phonons is negligible, due not only to charge cancelation 
but mostly to large LO phonon frequencies, compared to 
the timescales involved. 

Initially, the state of the system is described by the 
density matrix g = go <8) Q\{T), where go is the state 
with no exciton and g\(T) is the thermal equilibrium of 
lattice modes at the temperature T. The dynamics of 
the system is calculated by the perturbation expansion 
of the evolution operator in Born approximation. The 
lattice degrees of freedom are then traced out leading 
to the reduced density matrix for the carrier subsystem, 
which contains all the information accessible by optical 
methods. The details of the procedure may be found 
in2£. As a result, the reduced density matrix after time 
t may be written as g(t) = U(t)(g + gi)U^(t), where 
U(t) is the driven evolution for the unperturbed carrier 
system and g\ describes the perturbation to the carrier 



dynamics coming from phonons. The matrix elements of 
the latter, in the basis of empty dot (|0)) and one exciton 
(|1)) states, may be written as 

(0| ft |0) = -(l|^|l) = - f° du^SooM, (la) 



|f?i|i> = <iNo>* 



dto ^MsoiM.(lb) 



An additional component describing a unitary correction, 
e.g. light coupling renormalization and energy shifts, 
may be canceled by an appropriate modification of Hx 
(these effects may lead e.g. to intensity dependence of 
the observed Rabi frequencj^) . The function R(w) is 
the spectral density of the reservoir, fully characterizing 
the lattice properties at this order of perturbation treat- 
ment, defined as 



k 



l/fcl 2 [(n k + l)6(w - uj k ) + n k S(uj + uj k )] 



It depends on the material parameters and QD size via 
the coupling constants f k (see e.g^i for explicit formu- 
las). The functions Soo(w), S'oi(w) contain the complete 
information on the optically controlled carrier dynamics. 
If we assume that final state is measured after a time long 
compared to the phonon oscillation periods, they may be 
written as 

Soo(w) = \ [sin 2 a+|X s H| 2 ] , 

S*oiH = l[sm2a + 2ReK s (cj)K*(uj)], 
o 



where 



K,(u) = / dre luJT — sin2F(r), 

Js ",T 

K c {uj) = / dTe lu]T — cos2F(t), 
Js dr 



(2a) 
(2b) 



where s,i are the initial and final time of the evolution, 
F(t) = J dr/(r) is the rotation (on the Bloch sphere) 
performed on the exciton state up to time t and a is 
the total angle of rotation. For a given a, the functions 
(|2al b^l actually depend only on lut p , where r p is the pulse 
duration. 

In a typical experimenU2iiSii& one measures the aver- 
age dot occupation, (n) — (l|g(oo)|l), after a pulse of 
fixed length but variable amplitude. The simulation of 
such an experiment within the presented model is shown 
in Fig. njfor Gaussian pulses [r p is the full width at half 
maximum of the pulse envelope f(t)]. On the grounds 
of earlier numerical calculations^, we have assumed the 
exciton ground state to be of nearly product character, 
with Gaussian electron and hole wavefunctions with in- 
plane localization widths of 4.9 nm and 4.0 nm, respec- 
tively and with growth-direction width of 1 nm (this cor- 
responds to a dot of roughly 20 nm diameter and 4 nm 
height). The oscillations are almost perfect for very short 
pulses (~ 1 ps), then loose their quality for longer pulse 
durations (~ 10 ps). Although this might be expected 
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FIG. 1: (a-c) Pulse-area-dependent Rabi oscillations for var- 
ious pulse durations r p as shown in the figure, for T = 10 
K (a is the rotation angle on the Bloch sphere). Dotted line 
shows unperturbed oscillations, (d) The final QD occupation 
after two tt/2 pulses separated by time interval At for pulse 
durations as shown. 
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FIG. 2: Phonon spectral density R{u>)/uj 2 (solid) for deforma- 
tion potential coupling to LA phonons at T = 10 K and the 
nonlinear pulse spectrum Sqo{uj) (dashed and dotted lines) for 
t p = 1 ps, for rotation angles a as shown. 



from any simple decoherence model, the striking feature 
is that the effect dramatically grows for higher oscilla- 
tions, despite the fact that the whole process has exactly 
constant duration. Even more surprising is the improve- 
ment of the quality of oscillations for long pulses (~ 50 
nm) where, in addition, the first oscillation is nearly per- 
fect. 

The formulas (jlalb) quantify the idea of resonance be- 
tween the induced dynamics and lattice modes: in Fig. [21 
the phonon spectral density is compared to the non-linear 
frequency characteristics of the optically controlled exci- 
ton dynamics for r p = 1 ps (the characteristics for other 
durations is easily obtained by scaling). According to 
ifTalb) , the overlap of these spectral characteristics with 
the phonon spectral density gives the perturbation of the 
coherent carrier dynamics. 

For growing number of rotations n, the nonlinear pulse 
spectrum Sqq(lj) develops a series of maxima of increas- 
ing strength (Fig. |5Ji,b,c). The position of the last and 
highest maximum corresponds approximately to 27rn/r p , 
in accordance with the semiclassical resonance concept. 
However, spectral components are also present at all the 



frequencies 27rn'/r p , n' < n, which is due to the turn- 
ing on/off of the pulse. It is interesting to note that for 
high n, the low-frequency part of Sqo(lj) does not evolve 
with n (Fig. It is now clear that there are two 

ways of minimizing the overlap: either the pulse must be 
so short that all the maxima of Sqo(ui) are pushed to the 
right into the exponentially vanishing tail of the reservoir 
spectral density R(uj), or the pulse must be very long, to 
"squeeze" the spectral function near u> = and thus re- 
duce its area. In the latter case, the maxima developing 
with growing number of oscillations will eventually over- 
lap with R(oj) destroying the coherence. 

Although it might seem that speeding up the pro- 
cess is the preferred solution, it is clear that this works 
only because no high frequency features are included 
into the present model. In reality, speeding up the dy- 
namics is limited e.g. by the presence of excited states 
and non-adiabatically enhanced LO phonon coupling 28 . 
Moreover, it turns out that the resulting dynamics, even 
within this model, is actually not fully coherent. It has 
been shownS . that superposition of states created by an 
ultra-short tt/2 pulse becomes corrupted, preventing a 
second pulse (after some delay time At) from generat- 
ing the final state of (n) = 1 with unit efficiency. In 
order to prove the fully coherent character of carrier dy- 
namics it is necessary to demonstrate the stability of the 
intermediate state in a two-pulse experiment. The sim- 
ulations of such an experiment are shown in Fig. ^i. A 
short pulse (r p = 0.1 ps) creates a superposition of bare 
states (surrounded by non-distorted lattice) which then 
decohere due to dressing processe d 19 ^ 7 . As a result, the 
exciton cannot be created by the second pulse with unit 
probability^ .. For a longer pulse (r p = 1.0 ps), the lattice 
partly manages to follow the evolution of charge distri- 
bution during the optical operation and the destructive 
effect is smaller. Finally, if the carrier dynamics is slow 
compared to the lattice response times (r p ~ 10.0 ps), 
the lattice distortion follows adiabatically the changes in 
the charge distribution and the superposition created by 
the first pulse is an eigenstate of the interacting carrier- 
lattice system, hence does not undergo any decoherence 
and the final effect is the same for any delay time (its 
quality limited by decoherence effects during pulsing). In 
fact, splitting the 7r-pulse into two corresponds to slow- 
ing down the carrier dynamics which, in the absence of 
decoherence during delay time, improves the quality of 
the final state, as seen in Fig^i. 

The above analysis shows that damping of pulse-area- 
dependent Rabi oscillations due to interaction with lat- 
tice modes is a fundamental effect of non-Markovian 
character: it is due to a semiclassical resonance between 
optically induced confined charge dynamics and lattice 
modes. The destructive effect may be minimized both 
by speeding up and slowing down the dynamics. How- 
ever, in the former case, the system passes through un- 
stable (decohering) states. Moreover, fast operation on 
a real system induces many undesirable effects: transi- 
tions to higher states, bi-exciton generation or resonant 
LO phonon dynamics. On the other hand, for slow oper- 
ation, the number of "good" oscillations is limited. Thus, 
it is impossible to perform an arbitrary number of fully 



4 



coherent Rabi oscillations on an exciton confined in a 
quantum dot. 

By increasing the pulse duration as r p ~ a 2 , the 
phonon effect on the exciton dynamics may be kept con- 
stant. However, in this case the achievable number of 
oscillations is strongly restricted by the exciton lifetime 
and other (thermally activated) processes. Eliminating 
the radiative losses e.g. by using stimulated Raman adi- 
abatic passage instead of a simple optical excitation 29 
seems to be a promising direction from this point of view. 

The model presented above accounts for the decrease 
of the quality of Rabi oscillations in the short duration 
range observed in the experiment. It predicts, however, 
that this trend is reversed for longer pulse durations. The 
quantitative value of 96% for the first maximum of the os- 
cillations with t p = 1 ps agrees very well with the experi- 
mental result although the following extrema are much 
worse in reality than predicted here. This suggests an in- 
creased lattice response at higher frequencies which may 
be due to more complicated wavefunction geometr^S or 
electric-field induced charge separation leading to strong 
piezoelectric effectsi 9 -. Nevertheless, as far as it may be 



inferred from the experiment, the decrease of the oscil- 
lation quality seems to saturate after one full Rabi rota- 
tion, as predicted by the model calculations. 

In conclusion, we have studied the carrier-lattice dy- 
namics for optically induced Rabi oscillations of exciton 
occupation in a quantum dot. We have shown that the 
lattice response is resonantly driven by a combination 
of the (linear) pulse spectrum and the Rabi frequency. 
This semiclassical interpretation allows one to qualita- 
tively predict the damping effect based on the general 
knowledge of the spectral properties of the lattice modes 
and leads to the conclusion that the resonant lattice re- 
sponse precludes performing a large number of fully co- 
herent oscillations in this system. 
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